ON THE PROPERTIES OF THE COMBINATORIAL RICCI 
FLOW FOR SURFACES 



EMIL SAUCAN 

Abstract. We investigate the properties of the combinatorial Ricci 
flow for surfaces, both forward and backward - existence, uniqueness 
and singularities formation. We show that the positive results that exist 
for the smooth Ricci flow also hold for the combinatorial one and that, 
moreover, the same results hold for a more general, metric notion of 
curvature. Furthermore, using the metric curvature approach, we show 
the existence of the Ricci flow for polyhedral manifolds of piecewise 
constant curvature. We also study the problem of the realizability of 
the said flow in R 3 . 



1. Introduction and first results 

We consider the combinatorial Ricci flow introduced by Chow and Luo 
|13j . but we approach it not directly from its definition, but rather via the 
"classical" Ricci flow (i.e. the Ricci flow for smooth surfaces). While this 
approach is certainly less elegant, it allows us to obtain, by making appeal 
to the power of the classical theory (see [24] , [47] , [48] ) , a number of results 
that appear not to have been proved previously: existence of the reverse 
flow, uniqueness, singularities formation and the issue of embeddability (re- 
alizability) in M 3 . Such results are not less important from the viewpoint of 
the applications of the combinatorial flow - see, e.g. [23], [30] . 

We consider first compact surfaces without boundary and we concentrate 
mainly on polyhedral objects, since these are those arise more naturally 
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{S^} mg jj of smooth surfaces, (homeomorphic to Sp ol ), such that it assures 
good approximation of S P1 in the 



and also that are of most interest in applications. We shall also show that 
by considering a more general notion of metric curvature, we easily obtain 
an extension of these results to a larger class of geometric objects (and, in 
particular, to CW complexes, not necessarily reg ular)§ 

We start with the following basic fact: Given a polyhedral surface Sp ol , 
we can find a smooth surface S% arbitrarily close to it, both in the Hausdorff 
metric and as far as Gaussian curvature is concerned. More precisely, we 
have the following result due to Brehm and Kiihnel [8] 

Proposition 1.1 (Brehm and Kiihnel [5], Proposition 1). Let Sp ol be a 
compact polyhedral surface without boundary. Then there exists a sequence 

>PoV 

^Pol 

(1) metric sense, that is 

(a) S^ n = Sp ol outside the ^-neighbourhood of the 1 -skeleton of 
c2 

°Pol> 

(b) The sequence {S^ n } m& fq converges to Sp ol in the Hausdorff met- 
ric; 

(2) curvature sense, more precisely the (combinatorial) curvatures of S^ n 
weakly converge to Sp ol in the sense of measures. Here the curvature 
measure in the smooth case is the area measure weighted by Gauss 
curvature (considered as a function) and in the polyhedral case it 
is the Dirac measure at the vertices weighted by) the combinatorial 
curvature. 

Recall that combinatorial Gauss curvature for polyhedral surfaces is de- 
fined by the angular defect, that is: 

rrip 

(1.1) K(p) = 2Tr-J2ai(p) 

i=i 

where a\, . . . , a mp are the (interior) face angles adjacent to the vertex Vi. 
The combinatorial Ricci flow is then given as 

(1.2) ±± = -2K iri , 



lr The class of CW complexes for which this generalization holds is larger than the one 
considered in [13], since we do not have to restrict ourselves to complexes such that each 
to cell has 3 vertices and each vertex has degree 2. 
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in straight forward analogy with the classical flow 

(1.3) ^P- = -2K(t) 9ij (t). 

(In using proper, rather than partial derivatives notation, we follow [13j.) 
We also consider the close relative of (jl.3p . and its combinatorial version, 
namely in the normalized flow 

(1-4) ^ = -2(K - K(t)) 9lJ (t) , 

where K is the average sectional (Gauss) curvature of the initial surface Sq, 
K = f So K(t)dA/f So dA. 

Here, given the (Euclidean) triangulation with edges e^, of lengths lij, 
the numbers rj satisfy the following relation: 



(1.5) Uj = y r 4 2 + r| + 2r i r j cos($(e ij ) , 

where the r^'s represent the radii of the circle packing defined by the edge 
lengths above and by the intersection angles = Qifiij) between the 

adjacent circles of the packing. 

Note that, if starting with a polyhedral surface <Sp z, in order that (|1.3j) . 
Q1.4)l should hold (and, indeed, make sense as a flow), it is important to 
approximate Sp Q i both in the metric sense and also that curvature is well 
approximated, hence the interest in - and the need for - Proposition 11.11 

A proof of Proposition 11.11 is given in [8] (where a stronger result is also 
proved). In a nutshell, it proceeds as follows: Excise small enough (dis- 
joint) disk neighbourhoods of each of the vertices. The exterior of the union 
of such neighbourhoods is clearly approximable by a smooth surface com- 
posed of flat regions (the interior of the faces) and zero curvature cylinders 
(around the edges). Each of the neighbourhoods can be approximated by a 
smooth surface S(p), such that dS(p) is a geodesic m p -gon, such that the 
sum of its exterior angles is m p Tr — ^ Oj (p) . Obviously, as noted by Brehm 
and Ktihnel, the proof above extends to surfaces only locally embedded in 
]R 3 . (See also Remark 13.61 below regarding the existence of such local em- 
beddings.) Let us also note the fact that - a fortiori, given the freedom in 
the construction, provided by the supplementary dimensions - it also holds 
for surfaces embedded in M. N , for some N > 3% This is important if an 



2 A polyhedral model of the Boy surface [7] provides a good illustrating example 
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abstract surface is given and when one must start by considering first an 
isometric embedding in M. N . (See |41j . [42] for a discussion on the feasibility 
and practicability of such embeddings for PL, as well a more general types 
of surfaces, both in the applied and purely mathematical contexts.) 

However, we do not content with the previous demonstration, but we also 
provide a different proof, that also applies to a larger class of surfaces. More- 
over, it captures better the meaning of the notion of Hausdorff convergence 
and its interplay with curvature. 

Our approach is less elementary, but in a sense more natural for ge- 
ometers and topologists, as well as for image processing applications, and 
that uses only standard analytic tools: Instead of building the smooth sur- 
faces from a set of "standard elements" (cylinders, etc.), we consider in- 
stead smoothings (see, e.g. [33])§| Since, by [33], Theorem 4.8, such 
smoothings are 5- approximations, and therefore, for 5 small enough, also 
a-approximations of the given piecewise-linear surface Sp ol , they approxi- 
mate arbitrarily well both distances and angles on Sp ol . (Not to encumber 
the presentation with too many details regarding tools of differential topol- 
ogy, we have concentrated the relevant definitions, and results in an appen- 
dix.) Therefore angles, hence defects, are arbitrarily well approximated as 
well. While Munkres' results concern PL manifolds, they can be applied to 
polyhedral ones as well, because, by definition, polyhedral manifolds have 
simplicial subdivisions (and furthermore, such that all vertex links are com- 
binatorial manifolds) - see, e.g. jS], p. 346. The fact that 2-dimensional 
manifolds with a CW complex structure are also smoothable follows from 
the fact that any manifold of dimension < 3 admits a PL structure (see, 
e.g. [16]) and that, furthermore, this structure admits a unique smoothing 
(see, e.g. [33] )• In consequence, Gauss curvature of the smooth surface 
approximates arbitrarily well combinatorial curvature, at the essential com- 
mon points (i.e. the vertices of the given polyhedral surface). We note that, 
since, by [23], Corollary 5.2 (see also p3], Proposition 5.4) the Ricci flow is 
conformal, it follows that combinatorial curvature approximates arbitrarily 

'•A similar approach of approximating discrete structures by smooth ones is adopted 
also in theoretical physics [19], [20], the paradigm therein being that the structure of 
space-time at the smallest scales is, in fact, discrete and that classical models are smooth 
approximations of these structures. 
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well the curvature of the evolved surfaces, at any time t (see Section 2 below 
for details). 

Note that our proof renders in fact a somewhat stronger result than that 
of [8], since no embedding in M 3 is apriorily assumed, just in some H N ; 
however, as we have already seen, this represents only a slight improvement. 
More importantly, no change in the geometry of the 1-skeleton is made, not 
even in the neighbourhoods of the vertices. 

Moreover, it follows that metric quadruples (see definition below) on Sp a i 
are also arbitrarily well approximated (including their angles) by the corre- 
sponding metric quadruples) on S m . But, by [50] (see also |BM70| . Theorems 
11.2 and 11.3), the Wald metric curvature (see below) of S m , at a point p, 
K\y(p) equals the classical (Gauss) curvature K{p). Hence the Gauss cur- 
vature of the smooth surfaces S m approximates arbitrarily well the metric 
one of Spl (and, as in [8], the smooth surfaces differ from polyhedral one 
only on (say) the ^-neighbourhood of the 1-skeleton of Sp Q i. This state- 
ment can be made even more precise, by assuring that the convergence is 
in the Hausdorff metric. This follows from results of Gromov [2TJ - see |43j 
for details. That such curvatures converge not only punctually, but as mea- 
sures as well, i.e. that the so called CCP(i^) property of [8], Proposition 10 
also holds, follows, as a particular case, from [llj . Theorem 5.1, using the 
fact that polyhedral manifolds represent secant approximations of their own 
smoothings. 

Here metric quadruples are defined as follows: 

Definition 1.2. Let (M, d) be a metric space, and let Q = {p\, ...,pn} C M, 
together with the mutual distances: dij = dji = d(pi,pj); 1 < i,j < 4. 
The set Q together with the set of distances {djj}i<jj<4 is called a metric 
quadruple. 

Remark 1.3. Metric quadruples can be defined in a slightly more abstract 
manner, without the aid of the ambient space: a metric quadruple being de- 
fined, in this approach, as a 4 point metric space; i.e. Q = ({pi, ...,^4}, {dij}) , 
where the distances dij verify the axioms for a metric. 



this corresponds, essentially, to Condition (2) in Theorem 11.11 above 
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Before we can define the notion of embedding curvature, we have first 
to introduce some notation: Let S K denote the complete, simply connected 
surface of constant Gauss curvature k, i.e. S K = M 2 , if k = 0; S K = S 2 ^ , 
if k > 0; and S K = H 2 ^ - , if k < 0. Here S K = S 2 ^ denotes the sphere 
of radius R = and S K = H 2 // — - stands for the hyperbolic plane of 

curvature represented by the Poincare model of the plane disk of 

radius R = 1/ \J — k . 

Definition 1.4. The embedding curvature k(Q) of the metric quadruple Q 
is defined to be the curvature k of the gauge surface S K into which Q can 
be isometrically embedded. 

We are now able to bring the definition of Wald curvature [50] (or rather 
of its modification due to Berestovskii [1]): 

Definition 1.5. Let (X, d) be a metric space. An open set U C X is called 
a region of curvature > k iff any metric quadruple can be isometrically 
embedded in S m , for some m > k. A metric space (X, d) is said to have 
Wald- Berestovskii curvature > K iff for any i £ I is contained in a region 
of curvature > k. 

Remark 1.6. While the second part of the definition above is not needed 
in the remainder of the paper, we bring it for completeness (and for its 
importance elsewere - see, e.g. [H], [4"2]). 

Remark 1.7. Note that we can consider the Wald-Berestovskii curvature at 
an accumulation point of a metric space, hence on a smooth surface, by 
considering limits of the curvatures of (nondegenerate) regions of diameter 
converging to 0. 

Remark 1.8. The Wald-Berestovskii of a metric can actually be computed, 
using the following formula for the embedding curvature of a metric quadru- 









k{Q) = { 



K, K < 
K, K > 
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where dij = d(xi,Xj), 1 < i,j < 4, and D(Q) denotes the so called Cayley- 
Menger determinant: 



(1.7) 



Z?(xi,x 2 ,x 3 ,x 4 ) 
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However, it should be noted that, as far as the actual computation of k(Q) 
using Formula (|1.6|) is concerned, the equations involved are - apart from 
the Euclidean case - transcendental, are not solvable, in general, using ele- 
mentary methods. Moreover, when solving them by with the assistance of 
computer assisted methods, they display certain numerical instability. For 
a more detailed discussion and some first numerical results, see 1391 . 



Evidently, in the context of polyhedral surfaces, the natural choice for the 
open set U required in Definition 11.51 is the closed star of a given vertex v, 
that is, the set {e v j}j of edges incident to v. Therefore, for such surfaces, 
the set of metric quadruples containing the vertex v is finite. 

The combinatorial (defect) and metric (Wald-Berestovskii) notions of cur- 
vature are more closely related than just by having as limit, when the mesh of 
the polyhedral approximation of a smooth surface tends to zero, the Gauss 
curvature of the said (smooth) surface. Indeed, Wald-Berestovskii curva- 
ture can be characterized in terms of angles' sum at a vertex (hence defect). 
Before formally stating this result, we first need to introduce some further 
notation: 

Given three points X{ , Xj , X\ 111 3j metric space (X, d), we denote by a«(xj, Xj,x t ) 
£ [0, 7r], the angle Z(xjXjX/j^of the model triangle in S™. Let Q = {x±, X2, X3, X4} 
be a metric quadruple. We introduce the following quantity associated with 
Q: 

(1.8) V K (^Cj) — ct K (xj , Xj , x{) + cv K (xj , Xj , x m ) + cx K (xj , X[ , x m ) 

where Xj , Xj , x; , x m € Q are distinct, and k is any number. 

We can now bring the promised characterization of Wald-Berestovskii in 
terms of angle sum: 



■"that is of apex Xi 
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Proposition 1.9 (|36j. Theorem 23). Let (X,d) be a metric space and let 
U S X be an open set. U is a region of curvature > k iff V K (x) < 2tt, for 
any metric quadruple {x,y,z,t} C U. 

The result above shows that, in fact, the metric approach to curvature 
is equivalent to the combinatorial (angle-based) one, as far as polyhedral 
surfaces (in ]R 3 ) are concerned. (In fact, the metric approach is more general, 
since it can be applied to a very large class of metric spaces. Also, we 
should again emphasize that, as far as approximations of smooth surfaces in 
M 3 are concerned, both approaches render, in the limit, the classical Gauss 
curvature.) 

Note that we just gave a positive answer to the question - unposed so 
far, to the best of our knowledge - whether the metric curvature version of 
Brehm and Hiihnels's basic result also holds, namely we have proved: 

Proposition 1.10. Let Sp ol be a compact polyhedral surface without bound- 
ary. Then there exists a sequence {S^lmeN of smooth surfaces, (homeomor- 
phic to Sp i), such that 

(1) (a) = Sp ol outside the ^-neighbourhood of the 1 -skeleton of 

o2 

°Pol> 

(b) The sequence {5^} m£ N converges to S Pol in the Hausdorff met- 
ric; 

(2) K(S Pol ) — > Kyy(S Pol ), where the convergence is in the weak sense. 

Remark 1.11. The converse implication - namely that Gaussian curvature 
K(Y>) of a smooth surface S may be approximated arbitrarily well by the 
Wald curvatures Kw(Ep i im ) of a sequence of approximating polyhedral 
surfaces £p z, m - is, as we have already mentioned above, quite classical. 
(For other approaches to curvatures convergence, see, amongst the extensive 
literature dedicated to the subject, [H] and [6], [E], for the theoretical and 
applicative viewpoints, respectively.) 

We should also stress again the properties of the Gromov-Hausdorf conver- 
gence of finite e-neij^] in any sequence of approximating surfaces (polyhe- 
dral or smooth) of a given surface S 2 (again, smooth or not). In particular, 

"just for the record, recall that, given a metric space (X, d), a A C X is called an e-net 
iff d(x, A) < e, for any x G X. 
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by considering e-nets on surfaces, one automatically ensures (see [21], |10j ) 
any intrinsic geometric property of an approximating to the respective geo- 
metric property of the limiting geometric object. Most important for us, 
this holds for the intrinsic metric and thence for the metric curvatures - see 
[39J , @3j for a more detailed discussion and some numerical experiments. 

Also, one can consider simultaneously the combinatorial/metric Ricci flow 
on a polyhedral surface Sp ol , as well its classical counterpart on its smooth- 
ing S 2 . Since, as we have already noted above, the metric and the clas- 
sical curvatures, Kw(S Pol ) and K{S^ n ), respectively, are arbitrarily close 
to each other, and since the equations the two respective flows (combina- 
torial/metric and classical/smooth) contain the same curvature term, the 
ensuing metrics at each time during the flow will coincide on the common 
set, i.e the 1-skeleton of the polyhedral manifold and in the exterior of an 
arbitrarily small neighbourhood of it. Therefore, the limit surfaces for both 
flows - Sq Pol and Sg, respectively - will be isometric on the said set but 
perhaps only arbitrarily close to being isometric in the considered neighbour- 
hood!]] Moreover, by considering dense enough e-nets (of arbitrarily small 
mesh), the intrinsic metric of a polyhedral approximation gp Q i of a smooth 
manifold, and the (smooth) metric g of the later will be arbitrarily close to 
each other. (If one does not consider a limiting process, then the distortion 
of g by gp i can be computed, at each time "t" during the flow, using such 
formulae as (|2.2p and (|2.6p below, in conjunction with the computations in, 
e.g., [35], Lemma 3.19o This answers to a question posed to us by D. X. 
Gu. 

Before concluding this section, we should note that one can control the 
deformation both of the metric and of (combinatorial) curvature during the 
Ricci flow - for details see the proof of Theorem 12.41 below. 

2. Main results 

From the "good" , i.e. metric and curvature, approximations results above, 
it follows that one can study the properties of the combinatorial Ricci flow 

^One can ensure actual isometry by imposing a certain additional constraint on the so 
called "volume density" of the surface - for details, see (45 j . 

Q 

For manifolds with boundary, similar distortion estimates follow from the results in 
|40| . Of course, in this case, one still obtains isometry when restricting to the 1-skeleton. 
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via those of its smooth counterpart, by passing to a smoothing of the poly- 
hedral surface. The heavier machinery of metric curvature considered above 
pay off, in the sense that, by using it, the "duality" between the combina- 
torics of the packings (and angles) and the metric disappears: The flow is 
purely metric and, moreover, the curvature at each stage (that is, for every 
"i") is given - as in the classical context - in an intrinsic manner, i.e. solely 
in terms of the metric. 

A number of important properties now follow immediately. 

2.1. Existence and uniqueness (forwards and backwards). In partic- 
ular, the (local) existence and uniqueness of both the forward and backward 
(i.e. given by dgij(t)/dt = 2K(t)gij(t)) Ricci flows hold, on some maximal 
time interval [0,T];0 < T < oo (see, e.g. [J7], Theorems 5.2.1 and 5.2.2 
and the discussion following them). Beyond the theoretical importance, the 
existence and uniqueness of the backward flow allow us to find surfaces in 
the (conformal) class of a given circle packing (Euclidean or Hyperbolic). 
More importantly, the use of purely metric approach (based on the Wald 
curvature or any of other equivalent metric curvatures), rather than the com- 
binatorial (and metric) approach of [13], allows us to give a first, theoretical 
at this point, answer to Question 2, p. 123, of [13], namely whether there 
exists a Ricci flow defined on the space of all piecewise constant curvature 
metrics (obtained via the assignment of lengths to a given triangulation of 
2-manifold). Since, by Hamilton's results [23] (and those of Chow [T2], for 
the case of the sphere), the Ricci flow exists for all compact surfaces, it fol- 
lows from our arguments above that the fitting metric (combinatorial) flow 
exits for surfaces of piecewise constant curvature. In consequence, given a 
surface of piecewise constant curvature (e.g. a mesh with edge lengths sat- 
isfying the triangle inequality for each triangle), one can evolve it by the 
Ricci flow, either forward, as in works discussed above, to obtain, after the 
suitable area normalization, the polyhedral surface of constant curvature 
conformally equivalent to it; or backwards, to find the "primitive" family 
of surfaces (including the "original" surface) conformally equivalent to the 
given one. (Here, by "original" , we mean the surface obtained via the back- 
wards Ricci flow, at time T.) It is not necessarily true, however, that all 
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the surfaces obtained via the backwards flow are embedded (or, indeed, em- 
beddable) in M 3 - for details see Section 3 below. We can summarize the 
discussion above as 

Proposition 2.1. Let (S p gp a l) be a compact polyhedral 2-manifold with- 
out boundary, having bounded combinatorial (metric) curvature. Then there 
exists T > and a smooth family of polyhedral metrics g(t),t £ [0,T], such 
that 



(Here K(t) denotes the combinatorial (resp. Wald) curvature induced by the 
metric g(t).) 

Moreover, both the forwards and the backwards Ricci flows have the unique- 
ness of solutions property, that is, if g\(t), g2(t) are two Ricci flows on {Sp ol , 
such that there exists to G [0, T] such that gi(to) = 52(^0) > then gi(t) = g2(t), 
for all t £ [0,T]. 

In fact, the existence and uniqueness of the Ricci flow hold even if we do 
not restrict to compact surfaces, but we still require that the manifold is 
complete. Indeed, by applying the ideas employed in the proof above to a 
result of Shi [33] (see also [28])) we obtain the following 

Proposition 2.2. Let {S Pol ,gp i) be a complete polyhedral surface of metric 
(combinatorial) curvature bounded from above. Then there exists a (small) 
T as above, such that there exists a unique solution of 111. Sty for any t € [0, T] . 

Remark 2.3. Shi's result (hence the proposition above) does not necessarily 
hold for noncomplete surfaces - see |49j . 

Before proceeding further, we should stress that, at this point, the metric 
approach introduced above is purely theoretical and, while it allows for a 
number of (theoretical) results to be inferred from the classical theory, it 
lacks (at least for now) the simple algorithmic capability of the combinatorial 
one of Chow and Luo, and certainly of its subsequently development - see 



(2.1) 




23]). 
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2.2. Convergence rate. A further type of result, highly important both 
from the theoretical viewpoint and for computer-driven applications, is that 
of the convergence rate. For the combinatorial flow, it is shown in [13] that, 
in the case of background Euclidean (Theorem 1.1) or Hyperbolic (Theorem 
1.2) metric, the solution - if it exists - converges, without singularities, 
exponentially fast to a metric of constant curvature. Using the classical 
results of [23] and [12], we can do slightly better, since we already know 
that the solution exists and it is unique (see the subsection below for the 
nonformation of singularities). Moreover, we can control the convergence 
rate of the curvature: 

Theorem 2.4. Let (Sp ol , gp i) be a compact polyhedral 2-manifold without 
boundary. Then the normalized combinatorial ( metric ) Ricci flow converges 
to a surface of constant combinatorial (resp. metric) curvature. Moreover, 
the convergence (rate) is 

(1) exponential, if K < 0; x{S% i) < 0); 

(2) uniform; if K = 0; 

(3) exponential, if K > 0. 

Proof. As already noted, a unique solution for the Ricci flow exists for all 
< t < T and, again, these solutions are uniformly (conformally) equivalent. 
Indeed, by [21], Corollary 5.2 and [13], Proposition 5.15, there exists C = 
C(gp i), where gp a i is the smoothing of gp i, such that 

(2.2) -^gpoi <gt< Cgp i , 

where gt is the smoothing of gt , and the discussion above shows that the 
same holds for the polyhedral metrics. 

To estimate the convergence rate for the curvature, we make appeal to 
the following formulae (see, e.g. [13], Proposition 5.18): There exists C = 
C'(gpoi) > 0, (in fact, C = C'(gp i)), such that 

(1) If K < then 

(2.3) K - C'e Kt < K{t) <K + C'e Kt ; 

(2) If K = then 

(2.4) -_22_<* (t )< C < ; 
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(3) If K > then 

(2.5) - C'e Kt < K(t) <K + C'e Kt . 

To show that the metric also converges with exponential rate, one has to 
make appeal to a refinement of (12. 2p . namely that the constant C therein 
is, in fact, given by C = e 2KMax , where Kuax = max \K(t)\ , t G [0,T]. 
(This holds, in fact, for the case of the general Ricci flow, with Gaussian 
curvature K being replaced, of course, by the Ricci curvature Ric - see, e.g. 
|47j . Lemma 5.3.2.) In fact, given that the manifold under investigation is 
compact, hence of curvature bounded below and above, a stronger form of 
this improvement of (|2.2p can be given, and, moreover, one that is better fit- 
ted for the case of polyhedral manifolds (see, e.g. [31] Lemma 27.1, Remark 
27.5 and the following material): 

(2.6) e -« < ^4 < . 

dist (x,y) 

where Kuax is as above. 

By the approximation results above, namely Propositions 11.11 and 11.101 
the result follows for the combinatorial and Wald-Berestovskii curvatures, 
respectively. Alternatively, one can more directly infer the respective con- 
vergence rates from (|2.6p and (|2.2p as far as the metric is concerned, and 
for the Wald-Berestovskii curvatures from (|2.6p and (|2.3p _ (|2.5p . Similarly, 
the convergence for the combinatorial curvature follows again from (|2.6p in 
conjunction with (|1.8p and Proposition II. 9i □ 

Remark 2.5. Existence, uniqueness and, furthermore, convergence rate re- 
sults for polyhedral can be obtained, using the same techniques as before, 
for surfaces with (a finite number of) cusps and funnels, using quite recent 
results of Isenberg, Mazzeo and Sesum [28] (for finite area surfaces) and Al- 
bin, Aldana and Rochon pQ (for surfaces of infinite area). We do not bring 
here the technical details since they would bring us too far afield - for details 
and further related results, see 1291. 



2.3. Singularities formation. Another important aspect of any Ricci flow, 
be it smooth or discrete, is that of singularities formation. By [13], Theorem 
5.1, for compact surfaces of genus > 2, the combinatorial Ricci flow evolves 
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without singularities. However, for surfaces of low genus no such result ex- 
ists. Indeed, in the case of the Euclidean background metric, that is the 
one of greatest interest in graphics, singularities do appear [22] . Such singu- 
larities are always combinatorial in nature and amount to the fact that, at 
some t, the edges of at least one triangle do not satisfy the triangle inequality 
|22| . These singularities are removed in heuristic manner, using the graph- 
ics equivalent of e-moves (see, e.g. [E]). However, by [23], Theorem 1.1, 
the smooth Ricci flow exists at all times, i.e. no singularities form. By the 
considerations above, it follows that the metric Ricci flow also exists at all 
times without the formation of singularities. In fact, by a quite recent result 
of Topping [18], the same result holds even for unbounded (but complete) 
Riemannian 2-manifolds (M,g) with bounded curvature and satisfying a 
certain mild noncollapsing condition, namely that, there exists tq > 0, such 
that, for all x £ M, the following holds: 

(2.7) Vol 9 (B g (x,r )) > e > 0. 

(Here, as usual, B g (x,ro) denotes the open ball, in the metric g, of center x 
and radius tq.) 

Again, we can recap the discussion above as 

Proposition 2.6. Let (Sp j, gp i) be a complete polyhedral 2-manifold, with 
at most a finite number of hyperbolic cusps (punctures), having bounded 
combinatorial (metric) curvature and satisfying the noncollapsing condition 
\2. T\ j. Then there exists a unique Ricci flow that contracts the cusps. Fur- 
thermore, the curvature remains bounded at all times during the flow. 

Remark 2.7. Both the boundedness and the noncollapsing conditions evi- 
dently hold for surfaces that appear in graphics, hence the fitting result for 
the combinatorial (and metric) flow also hold for this type of application. 
It follows, that we can apply also in this context the cusp contracting result 
of 08]. One may, however, argue that manifolds with (hyperbolic) cusps do 
not appear in graphics, only compact manifolds (with or without boundary), 
but in fact many algorithms are modeled upon surfaces with punctures - 
see e.g. [23], [30] . 



3. Embeddability in R 3 
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In this section we mainly consider a problem regarding smooth surfaces, 
and we hope that by now, the connection with its version for polyhedral 
surfaces is clear. It should be noted, in this context that, by [33] . Theo- 
rem 8.8, any 5- approximation of an embedding is also an embedding, for 
small enough 5. Since, as we have already mentioned, smoothing represent 
5- approximations, the possibility of using results regarding smooth to infer 
results regarding polyhedral embeddings is proven. (The other direction - 
namely from smooth to PL and polyhedral manifolds - follows from the 
fact that the secant approximation (see Appendix) is a ^-approximation if 
the simplices of the PL approximation satisfy a certain nondegeneracy con- 
dition - see [33], Lemma 9.3.) We wish to stress here the importance of the 
embeddability in graphics and image processing. In the only fully imple- 
mented Ricci flow, that is the combinatorial flow [23], [30], the goal is, in 
fact, to produce, via the circle packing metric, a conformal mapping from 
the given surface to a surface of constant (Gauss) curvature. Since in the 
relevant cases (see [13]) the surface in question is a planar region (usually 
a subset of the unit disk), its embeddability (not necessarily isometric) is 
trivial. Moreover, in the above mentioned works, there is no interest (and 
indeed, no need) to consider the (isometric) embeddability of the surfaces 
Sf (see below) for an intermediate time t ^ 0, T. However, this aspect is 
very important if one considers the problem of the Ricci flow for surfaces of 
piecewise constant curvature; as well as in image processing - see [3]. 

Let Sq be a smooth surfaces of positive Gauss curvature, and let Sf denote 
the surface obtained at time t from Sq via the Ricci flow. For all omitted 
background material (proofs, further results, etc.) we refer to [25] . 

Proposition 3.1. Let Sq be the unit sphere § 2 , equipped with a smooth 
metric g, such that K(g) > 0. Then the surfaces S$ are (uniquely, up to a 
congruence) isometrically embeddable in M 3 , for any t>0. 

Proof. By a by now classical result of Nirenberg [M] and Pogorelov |37j 
(independently), Sq admits a smooth isometric embedding in M 3 . 
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The metric gt being conformal to go, for any t (see [21]), we can write it as 
gt = e 2ip go, for some smooth function ip% Then, by elementary computations 
(see, e.g. [25]), we obtain that the Gauss curvature of this metric is 

K t = K{g t ) = e-^{K go - tA go <p) = te^^K, + (1 - t)e~ 2t * > . 

Therefore, again by Nirenberg and Pogorelov's theorem, the surfaces S 2 are 
isometrically embeddable in M 3 , for any t > 0. □ 

We can, in fact, do somewhat better: 

Corollary 3.2. Let S 2 be a smooth surface. If x{Sq) > 0, then there exists 
some to > 0, such that the surfaces S 2 are isometrically embeddable in K 3 , 
for any t > to. 

Proof. By the continuity of the Gauss curvature during the Ricci flow, it 
follows that, as some time to, the \Kt \ < Kq > 0. Applying again the 
arguments in the proof above, the corollary follows. □ 

Remark 3.3. In this context it is impossible not to mention Alexandrov's 
results [2] regarding convex surfaces in M 3 : (A) Any convex surface, en- 
dowed with its intrinsic metric, is a manifold of nonnegative curvature; and, 
in essentially the opposite direction, (B) Any complete metric of positive 
curvature (nonnegative) on the sphere (the plane) represents the metric of 
a (not necessarily smooth) convex surface. 

However, we should underline that requiring that a certain polyhedral 
sphere actually has (strictly) positive curvature at all its vertices it is quite 
a strong condition; indeed, it is a well known fact in graphics (see, e.g. 
|32| ) that even the most standard polygonal approximations of the sphere, 
exhibit, even at and high resolution, saddle points at certain vertices. 

In contrast with this positive result, the for (complete) surfaces uni- 
formized by the Hyperbolic plane we have the following negative result: 

Proposition 3.4. Let (S 2 ,go) be a smooth surface. If \{S 2 ) < 0, then 
there exists some to > 0, such that the surfaces S 2 are not isometrically 
embeddable in ]R 3 , for any t > to- 

^This represents, in fact, the first, basic step in the proof of the Niremberg-Pogorelov 
theorem. 
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Proof. Since at time T, the surface undergoing the flow has constant negative 
Gauss curvature, it is not smoothly (C 4 ) embeddable in M 3 , by a classical 
theorem of Hilbert [27]. By the continuity of the Gauss curvature during 
the Ricci flow, it follows that, as some time to, Kt < Kq < 0. Therefore, by 
a result of Efimov [T7J, it follows that 5( admits no smooth (C 2 ) isometric 
immersion (hence embedding) in M. 3 . □ 

Remark 3.5. Efimov \18\ also proved that even if only the gradient of the 
Gaussian curvature is bounded (by some specific constant - see [H], [25]). 
there exists no smooth (C 3 ) isometric immersion in M 3 , and no C 2 isometric 
immersion exists if the suprema of K and of its gradient are < oo. On the 
other hand, Hong [26] showed that smooth isometric immersions in M 3 exist 
if the decay rate at infinity of K is slower that the inverse square of geodesic 
distance. 

Obviously, the above mentioned results render appropriate versions of 
Proposition 13.41 However, their applicabilty is of far lesser interest in the 
context of the Ricci flow for polyhedral surfaces, so we do not formulate 
them explicitly. 

Remark 3.6. In fact, the the situation is far worse, so to speak, than even 
Hilbert's and Efimov's theorems might suggest. Indeed, no general exis- 
tence/nonexistence result is available, even if one restricts himself to local 
isometric embedding. Without going in too many details (since this would 
bring us out of our scope), it is known that such embeddings are possible if 
K does not vanish; or when K{p) = and dK(p) ^ or if dK(p) > in 
some neighbourhood of a point p; and again when K{p) = 0, alK(p) = and 
HessA'(p) < 0. (For further details and bibliographical references see e.g. 
[25].) 

However, as in the global case, (for lower differentiability classes), no 
general results are even possible. This was first shown by Pogorelov [38] . 
who constructed a C 2 ' 1 metric on the unit disk B 2 = B 2 (0, 1) C R 2 , such that 
there exists no C 2 isometric imbedding in M 3 of B 2 (0, r), for any < r < 1. 

We mention in this context that a criterion for the local isometric embed- 
ding of polyhedral surfaces in M 3 , akin to the classical Gauss fundamental 
(compatibility) equation in the classical differential geometry of surfaces, 
was given in [41] . Namely, given a vertex v, with metric curvature Kw(v), 
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the following system of inequalities should hold: 

{maxAo(v) < 27r; 
a (v;vj,vi) < a (v;vj,v p ) + a (v; v h v p ), for all Vj,vi,v p ~ v; 
V K {v) < 2tt; 

Here 

(3.2) A) = maxVb; 

i 

"~" denotes incidence, i.e. the existence of a connecting edge e, = m/j 
and, of course, V K (v) = a K (v;vj,vi) + a K (v; Vj, v p ) + a K (v; v h v p ), where 
Vj,vi,v p ~ v, etc. 

Note that the first two inequalities represent the (extrinsic) embedding 
condition, while the third one represents the intrinsic curvature (of the PL 
manifold) at the vertex v. 

For details and a fitting global embedding criterion see [41j . 

We mentioned these facts because, beyond the "archival" interest (so to 
speak) in such results, they are relevant if locally embedding based applica- 
tions are envisioned (even though, as we have previously shown, the basic 
smoothing and approximation result holds in this case as well). 
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Appendix 

We include here the modicum of differential topology needed (mainly in 
the alternative proof of Proposition [TTT]). Our source for this material is [33j . 
We presume that the reader is familiar with the basic concepts (simplicial 
complexes, triangulations, etc.) however, as a background text, we warmly 
recommend Munkres' notes [33]. 

Definition 3.7. (1) Let / : K -> W 1 be a C r map, and let <5 : K -> 
be a continuous function. Then g : \K\ — > IR n is called a 5- 
approximation to f iff: 

(i) There exists a subdivision K' of K such that g £ C r (K' , M n ) ; 
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(ii) d eucl (f(x),g(x)) < S(x) , for any x G \K\ ; 

(iii) d eud (df a (x),dg a (x)) < 5(a) ■ d eud (x,a) , for any a G \K\ and for 
all x G Si(a,K'). 

(2) Let K' be a subdivision of K, U = U, and let / G C r (K,M n ), 5 G 
C r (i^',]R n ). g is called a 5 -approximation of / (on {/) iff conditions 
(ii) and (iii) above hold for any a G U. 

(Here St(a, K) denotes, as it standardly does, the star of a (in 
K),i.e. St(a,K) = {J a&a ^ K a. 

Recall that in the PL context the differential (of a map) is defined as 
follows: 

Definition 3.8. Let a be a simplex, and let / : a — > R n , / G C r . If a G a we 
define df a : a — > W 1 as follows: df a (x) = Df(a)-(x — a), where Df(a) denotes 
the Jacobian matrix Df(a) = (dfi/dx^i^ij^n, computed with respect to 
some orthogonal coordinate system contained in n(cr), where n(cr) is the 
hyperplane determined by a. The map df a : a — > W 1 does not depend upon 
the choice of this coordinate system. 

Moreover, df a \ ariT is well defined, for any a, r G St(a,K). Therefore the 
map df a : St(a, K) — > M n is well-defined and continuous, and it is called - 
analogous to the case of differentiable manifolds, the differential of /. 

Definition 3.9. Let K' be a subdivision of K and let / G C r (K,R n ), g G 
C r (K',M n ) be non-degenerate mappings (i.e. rank(f\ a ) = rank(g\ a ) = 

o 

dim a, for any a G K) and let U = U C \K\. The mapping g is called 
an a- approximation (of / on U) iff: 

(3.3) Z[dfa(x), dg a (x)) < a ; for any a G U, and any x G St(a, K'), a/i. 

As expected, a fine enough (5-approximation is also an a-approximation: 

Lemma 3.10 ([33], Lemma 8.7). Let K be a (finite) simplicial complex and 
let f : K —> M. n be a non- degenerate C r , 1 < r < 00 map. Then, for any 
a > 0, there exists 5 = 5(a) > such that any non- degenerate C r map 
g : K' — )■ M n , , which is a 5-approximation of f on some open set U, is also 
an a -approximation of f on U . (Here K' denotes, as above, a subdivision 
ofK.) 

We conclude the appendix with the following definition: 
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Definition 3.11. Let / £ C r {K) and let s be a simplex, s < a £ K. Then 
the linear map: L s : s — > W 1 , defined by L s (v) = f(v) where v is a vertex of 
s, is called the secant map induced by f. 
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